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Abstract 

The deep-water wave problem of any initially prescribed periodic stationary surface 
pressure has been discussed. Asymptotic analyses of both the steady and the unsteady parts 
of the resulting wave integral are given. Exact solutions for the steady state waves have 
been obtained for the cases of concentrated and arbitrary symmetrical pressure distributions of 
finite extent. Some particular Bpatial pressure configurations over elliptic regions together with 
some characteristic features of the wave motion are discussed. 

Introduction 

The three-dimensional problem of deep-water surface waves due to a 
stationary periodic surface pressure has been investigated by Lamb ( 1904 ), 
Sretenskii ( 1956 , 1957 ) and Wehausen and Laitone ( 1960 ). Of these, Lamb 
obtained the steady state integral for the wave elevation due to any symmetri¬ 
cal pressure and asymptotically evaluated this integral for the case of a con¬ 
centrated pressure point. Investigation of Sretenskii was mainly concerned 
with the steady waves from a doubly modulated pressure distribution over a 
rectangular region, while Wehausen and Laitone gave only the velocity 
potential integrals for the general initial value problem with an arbitrary law 
of time-dependence and also for the steady state problem with harmonic 
pressures. In the second case, a first order asymptotic expression for the 
potential integral and the average rate of w r ork done by the pressure were also 
given for symmetrical pressure distributions. 

In this paper, the deep-water wave problem of any initially prescribed 
stationary periodic surface pressure has been solved by a method different 
from that used by Wehausen and Laitone. This method is equally applicable 
to the case of arbitrary time-dependence of the initial surface pressure. The 
following new points of.interest are discussed. Firstly, methods of finding 
complete asymptotic expansions of both the steady and the unsteady parts of 
the resulting wave integral are given for different ranges of time and distance. 
The steady state wave integral has been exactly evaluated for concentrated and 
arbitrary symmetrical pressure distributions. An expression for the mean 
rate of work by an arbitrary pressure on the fluid is given. In all these cases, 
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particular spatial pressure configurations over elliptic regions, hitherto unin- 
vestigated, are studied. Some characteristic features of both steady and 
unsteady waves are discussed. 

Problem 


A semi-infinite ocean of a heavy homogeneous inviscid liquid is initially 
at rest. Its undisturbed horizontal free surface forms the #^-plane "while the 
s-axis is taken vertically upwards. Surface waves are excited and maintained 
in the ocean by the application of a surface pressure po(x, y, t) of the form 

p 0 (x i y > t) = F(z ) y)e**, t > 0, .(1) 

in which the arbitrary function F(x, y) is supposed to possess a double 
Fourier transform. If <j>(x , y, z; t) be the velocity potential of the irrota- 
tional motion so generated, the equation of continuity is 

V 2 <£ = 0, z < 0, * > 0. (2) 

Introducing a frictional force equal to y times the velocity, and assuming 
small oscillations in which squares of the velocities are negligible, we have, 
from Bernoulli's pressure equation, 

|+f/2-/4- ^ = F^t), 


where p denotes the pressure at [x, y, z) at time t, p the uniform density of the 
liquid, g the acceleration due to gravity and F^t) an arbitrary function of 

dd> 

time t. On absorbing F^t) in ^ , the surface condition (1) gives 

gpi= —F(x t y)e iwt + p^fi4> + ^ .(3) 


where £ denotes the surface elevation at the point (x, y. 0) at time t above the 
undisturbed surface. 

The kinematical surface condition gives, to the above degree of 
approximation, 



, *> 0 . 

o 


(4) 


Eliminating £ between (3) and (4), we get 




(5) 


The conditions at infinity are 

(f> and its first and second derivatives 0 

as 2 ->—oo, t > 0. .. .. .. (6) 


The initial conditions are 
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^(x 5 y, 0 ; 0) = 4> t (x, y, 0 ; 0) = 0.(7) 

The problem is to find a solution <f>{x, y, z; t) of (2) which satisfies the equa¬ 
tions (2), (5), (7) together with the condition (6). The surface elevation £ is 
then obtained by (3). 

Solution 

Assuming a double Fourier transform <J>(K l9 K 2 , z; t) for the velocity 
potential function y, z; t) namely 

fkKi, Kf>, z ; t) = ± J“ J 4(x, y, z ; t)e^ x+K ^dx dy, . . (8) 

we obtain from (2) the differential equation 

g-(jrf+A'^ = o. 

Therefore 

$ = A(K U K t , i)e‘(** +K *y l * .(9) 


in view of (6). 

The result of transforming (5) is 


d 2 J> d& d<l icv = 

w* +>l i + < l i = T m ” 


in which F(K l9 K 2 ) is the double Fourier transform of F(x y y). Using (9), 
we find that A satisfies the differential equation 


-jF(K lt K 2 )e iwt .( 10 ) 

e o2 = gK&ndK=(K 2 1 + Kl) i .(11) 

On putting A = v(K^. K 2 , t)e~ i,a and neglecting y 2 , we obtain from 


where 


(10) the equation 


™ +a 2 v = ^ e MF {Kl ,K 2 ). 


The real solution of this equation is 

v = v 0 {K u K 2 ) cos {at-a 0 (K lt K 2 ) } +— F(K lt K 2 ) f e (i 

P° Jo 

The initial conditions (7) make Vo — 0. 

Consequently, 

A = — F(K v K 2 )e i0>t f sin as 6 

P° J n 




sin o(tr-s)ds. 


■■ ( 12 ) 
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Inserting the last expression for A in (9), and inverting (8) by the Fourier 
double inversion theorem, we find that the velocity potential is 


<f>(x,y,z;t) = ~e ia “ f f - F(Ki, f e _(i ' i+< " )s 

j - x j _ x CT J 0 


X sin as ds dK x dK z 


(13) 


If we put /x = 0 in (13), the result agrees with that of Wehausen and Laitone 
(1960, 22.48) for the case of harmonic strengths of the initial pressure. 

Equation (3) now gives, after a little simplification in which /x 2 has been 
neglected, the expression for the surface elevation as 


<Jpi = 


-F(z, y)e i(ut + 


sJJl 


F(K lt Ko) e Kz_i(K i' t+K 2'' ) 


1 -}- ifuo 


-m2 


X sin at— sin at— (fx+icj) cos crjj dK x dK 2 

**JJ“ F(K„ K 2 ) «-<(*»*+V) 


ificu—w 2 ia)t 


( 7 2 + iflOJ — co 2 


dK x dK 2 . 


(14) 


The passage to the limit z -> 0— in the second term of (14) is easily justified for 
finitely extended pressure distributions F(x , y)e' wt . For a concentrated dis¬ 
turbance, this is justified by the results obtained. 

Finally, the solutions of the present problem are to be obtained by passing 
to the limit 0+ in (13) and (14). 

We can write 

£ = £i + £2 > • • * • * • * • (15) 

where gpl x and < 7 p£ 2 are respectively the second and the sum of the first and the 
third terms on the r.h.s. of (14). Clearly, Ci represents the transient part of 
£ which dies out (even in the limit p ->• 0 + , as wall be shown later) as t -> oo, 
while £ 2 represents the steady state oscillatory solution for the wave ele¬ 
vation £. The general solutions of the less general problems of (i) symmetrical 
and (ii) certain asymmetrical oscillatory pressure distributions over circular 
regions follow at once from (13) and (14). 


(i) 


Symmetrical pressure areas 
If F{x, y) = /(r), 


r = (z 2 +y 2 )*, then 


F(Ki, K t ) = f(K), K = (Kl+K:)\ 

where j(K) = j a/(a)/ 0 (i^a) da is the zero-order Hankel transform of the 
J 0 

function f(r) (Bochner and Chandrasekharan 1949, Th. 40). By the same 
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theorem, the velocity potential <j>(r, z; t) and the wave elevation £(r, t) in this 
case are found to be given by 

4>{r, z; t) ==— f — j(K)e Kz J 0 (Kr) f e~ <ili+iw)s sin as ds dK .. (16) 

P J 0 a j 0 


2 


Ke Kz f(K)J 0 (Kr) 
a 2 +ifjuo — (o 2 


gp((r, t) =* -“f(r)e iwt +ia)e z lim f 

z — > 0 — J 0 

X J^cr sin at— sin at— (M+uo)'cos <rtfj dK 


(17) 


Combining the first and the third terms on the r.h.s. of (17), we find the part 
£ in the steady state to be 



K 2 ](K)J 0 (Kr) dR 

aZ+ifiaj — cd 2 


of 


This is Lamb’s well-known steady state integral for the wave elevation in 
this problem (i) (Lamb 1904, eqn. 123). 


(ii) Asymmetrical pressures on circular regions 

If F(x , y) — fi(r) cos md+f 2 (r) sin md , 


where r, 6, z are cylindrical co-ordinates, and m is a positive integer, it can be 
easily shown that 

F(K lt K 2 ) = e imnl2 [J 1 (K) cos m>P+f 2 (K) sin m,p] 

where K, 0 are the plane polar co-ordinates in the (iCi, ijL 2 )-plane, and ji(K) 
and / 2 (-£) are the mth order Hankel transforms of the functions/i(r) and /. 2 (r) 
respectively. 

In this case, the velocity potential <j>(r 3 6, z; t) and £(r, 6\ t) are found 
to be given by 

4>(r, e, z ; t) = — e** f F 12 (K, 9) - e K ‘J m (Kr) f e' ( ^ + * u,s sin as ds dK (18) 
P J o ° Jo 

gpl(r , d;U) = — [/i(r) cos mO+f 2 (r) sin md]e i<ot 


+iwe 2 lim 

z-^0- 



KF l2 (K, 9)e Kz J m (Kr) 

a 2 + ifJLco — a > 2 


X [or sintrf— 


ifiw 

2a 


sin at—(fi+ico) cos at) dK 


+ (i/z u>—w 2 )e twi 



KF^K, 6)J n (Kr) 

cr 2 -j-ifx.w —to 2 


dK .. (19) 
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where Fn{K, 6) = f^K) cos mO+] 2 {K) sin mO, .. .. (20) 

In (16) to (19), the final expressions are to be obtained by passing to the 
limit M 0 + . 


Asymptotic analysis of l\. 

A method of finding the complete asymptotic expansion of the transient 
part £ 1 of £ (in the limit M -> 0 +) for a certain range of large values of r and t will 
now be given. We assume that the impressed pressure p Q (x , y, t) is zero 
outside a finite region D of the xy -plane. 

As defined in (15) 


gpt i = sr e 2 lim 

ZTT ar —0 — 


F(K lt K 2 )e Kz ^ K I x+K 2 y ) 

—o>2 


X I a sin ot - - sin cr^— (fi + ita) cos cr^J dK\ dK 2 

in which finally the limit 0+ is to be taken. 

Let 


/ = lim f f F(K U K 2 )e K ‘- i ( K i* +K z y '>. — 

Z-+0- J J 0(0 Z + l 


-\-i}XCO —6U 2 ) 


dKi dK 2 . .. (22) 


Replacing F{K U K 2 ) by the corresponding double integral and then inter¬ 
changing the orders of the two double integrations, we get 


=// 


F(ct, P)da. dfi .Ii. 


so that 


*'■ = ,lit fJJ^' 

-oo 

= lim f e K ‘J 0 (KR) 
o- J n 


l( ,-« )+W)] _^LZ^. dKidKz 


Kb j / tt t\\ K sin at j _ 

' J °< KB) «K=lF) dK 


where 


.R2 = (x—<t) 2 +(y—fi) 2 

rr, ft( 2 —t/MO 


Replacing J 0 (KR) in (24) by its integral representation 


J^KR) 


75 

_ 2 p 


cos (KB sin y) dy 
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and then interchanging the orders of y- and if-integrations, we get 

T>gh = (~ 2 dy lim f Ke Ks cos {KR sin y) — ^ dK. 

* J o *-*0- J o °A-n “"A ) 

Expanding sin at in powers of at, and interchanging the orders of summation 
and ii-integration, we obtain 


J 0 * = o v 


+TFT.!!?- 


A"+l 

"g_j£i &K cos {dili sin y) dK 


2 rfy 2 (- 1 )” 2 ^ ]i “ Z ”"'«** cos (KR sin y) rfA 


n«0 


y = 0 


rfy 

o n = 0 


(£)‘*■>(»*'<*)* 

S'- 1 >" (SfTTi 2 K " CL * <"-'+» I 


1-2 sin (gK'l 2 )* 


where 
12 
Since 


- 1 (f) 

-/!»*■■. 

1 £(2/+l, 2n-2r+l) |^ 2 " 2 ’(l-s) 2r &' 


( 20 ) 


(2n+1)! r(2r+l)r(2»-2r+l) r(2r+l)r(2n-2r+l) ’ 

we can apply Cauchy’s rule for the multiplication of power series, and obtain 


^ j 1 2 sin Tjj = 4 {— j .4 J ‘ cosec Y dv J cos r'(1 —s)L(2t 2 s 2 ) d$ 

where we define the functions L and M by 


(27) 


T , , •* , x Zn+1 2 Zn (2n )! 

I'M = J o (-!)•... ( 5 ?iTr 

() ~„to ( } (4w+3)! 


•• (28) 


and some dimensionless variables r etc. by 


g * 2 


gt 2 


t = 47 ’ T 1 T 2 = T! cosecy 


r' = tut, t 2 = (gK't 2 )i, r' = — 


•• ( 29 ) 
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where r 2 = (* 2 +/3 2 ) max , (a, /5) being a point in the domain of F(x, y). 

It is known that 

L(x) = {nx) 1 j^cos | C (|) + sin ~8 j.(30.1) 

= (ttz) 4 |cos |+ sin J e~^ 2x sin A 2 <Za| , x > 0 .. (30.2) 

and 

M(x) = [sin | (7 -cos |,S' (j)]..(30.3) 

(Ryshik and Gradstein 1957, 6*233 (2, 3), 6-236 (3)) where C(x) and S(x) 
denote Fresnel’s integrals defined by 


[C(x) 9 S(z)] 


= j_r 

( 2 »)‘ J „ 


[cos t y sin t] 


dt. 


By (30.2), 


the s-integral of (27) = J s cos t' x (1 —s)[eos (r/)+siii (V 2 )]d« 


—2t| f sin X 2 d\ f 

^ o J 0 


s cos rj (1—s)c ^ 2 ds 


(31) 


in which the orders of s- and A-integrals have been interchanged in the second 
term. 

In the asymptotics that follow, the range of large values of r and t is 
defined by the conditions: 

r'> 1, .(32) 


Applying Erdelyi’s theorem (Erdelyi 1955, 2.9) to the first term of (31) in 
which r 2 is predominantly large, we have now the first term of (31) 

= “ [cos|T 2 -(«+l)f|+sin{T 3 -(«+l)||] 



+ 1-JU 



COS Tj 



.. (33) 


where 


l(v) = — i cos {t'(1 — y/l— v)} 


.. (34) 
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Next, evaluating the ^-integral of the second term of (31), we get 
the second term of (31) 


= 2r* jV^si 


• (2Ar|+l)(4r 2 A 2 +r; 2 )-2r; 2 

sin A 4 •-:-7777~o-oA 


(4r^ +r y 

—2 t |[(/ 8 —2t' 2 / 4 ) cos Tj+ 4sin r'] .. 


( 35 ) 


where 


h, 1 


h/ 5 = r 

j 0 


(1, 1, A) sin A 2 


dX , taken in order. 


( 4 r 2 A 2 +T ; 2 )^’ 2 

Since r 2 ^> 1, the asymptotic expansion of the first term of (35) can be com¬ 
pletely found out by Watson’s Lemma (Watson 1944, § 8.3). Using some 
results on Stieltjes’ transforms (Erdelyi 1954. 14.2 (31, 33, 48) and certain 
properties of sine and cosine integrals, we can write 

1 r r cos (KR sin y) . r'l r*' 

h — -» - dK+in cos- \-tt sin — 

32t 2 [J o K-K' 4 T2 4tJ 

while 1 1 is easily evaluated from I 3 . 

Taking all these results for (35) together with (33), we find the s-integral 
of (27) 

= "(-?) 2 7^ [ cos j T 2- (w+1 )l( + dn { t 2-(»+ 1 ) if] 


n = 0 '2 


\ T i t 2 / Jo 


cos (KR sin y) 
K-K' 


dK+-~ le- iT i 0 os —. 
4 t „ 4 t ' 


.. (36) 

It is now clear from (27) that the complete asymptotic expansion of /, (and 
hence of / and f i) can be found out, the y-integrals involved being calculated 
as shown below. However, retaining only up to the second order terms, we 
get, from (27), 


+3rd order terms 

= it | —I I {^( T i cos fi A)— cosh A)} dA+3rd order terms 

\ 9 J t 1 J o 

- - ¥ (f)'? k (?) r > (?) (?) T -> (?) !+ m « d - 


(Erdelyi 1963, 7.14.2 (66)). 
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Since t, is large, we may use the asymptotic series for the Bessel and Neumann 
functions written above. Then 

w r r lK'\ k -ir' j 7T (2K'V/ . 3 \ , „ , , 

2 ^i — 77 J e VofA i?)J = — ^ smT i“ g ~ cosTij+Srd order terms. 

Expanding B~ l in tj in powers of r" 1 and retaining only those terms which 
contribute up to the second order terms, we have, from (23), 

I = J^K ( +n e~ il v '*' J F(a, p)J 0 (K'R) da dfi 

+3rd order terms 

where 3f m (z) ~ the imaginary part of 2 , 


• (37) 


F = F(w cos 6, w sin 6), w = 


2gt 2 
4r 2 


p i_, 1 52 

r ww ~ Sw 2 ’ V “ ~ 3,„2 + ,„ a«i + 


The other members of (21) being proportional to ^, we easily find 


Bw 2 ~w dw~w' 2 ' d$ 2 ’ 

BI d 2 1 
dt ’ Be¬ 


rn 


— ( W z e * ) +2w ft, ( Wtf * 


Tgr 


under the restrictions (32), 
where 


igt* 


& 


+3rd order terms (39) 


Wi = W+ ^ F ~ l K ’ Wi = w '- t Fm> 


(40) 


W being given by (38). 
The lowest order term is 


“ Wgr [ 2a,R ‘ ^ ( Fe ^)] ’• •' (41) 


under the restrictions (32). 

When F(x , y) = f(r) y F(w cos 6 , w sin 0) = /(w), 

where ]{k) is the zero-order Hankel transform of/(r); then 

w i «/(»)-» 

W z = W 1 - l ; ?(w) 
and £i is obtained from (39). 


( 42 ) 
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When * 

F{z, y) —/i(r) cos mQ+fe{r) sin mO , 

F(w cos 0 , w; sin 0) = e tm7T/ ' 2 cos m^+/ 2 (^) sin mO], 

fi{K) and / 2 (i^) being the mth order Hankel transforms of f^r) and f 2 (r) 
respectively. The corresponding expression for fj follows from (39) by (38) 
and (40). 


Decay of the unsteady state 
By (24) 


/•CO 

l x = lim G 
z-* 0 - J o 


G(Jc) sin (gKW)t dK 


G( R) « J 0 (K*It)eK 2 ‘ . 

The integral is a Fourier integral and the function G(K) satisfies all the 
conditions required for this integral to have an asymptotic expansion as t-> oo, 
r and o> remaining fixed. 

Therefore, we have (Erdelyi 1955, 2.8) 




Consequently, from (21), 

gPii = l £e~[-(2n+l)(2n+2) + (r+iw)(2n+l)t-^t^ 
• (2n+l)M-(2 re +2)], as M 0+, .. 


=i t 


where 


To find 


a„ = (-1 y2g-< n+ *> J J F(a, 0) jlim_ G {2n \ 0) J doc dp. 


lim 6r (0), we note that 

Z->0~ 


G (2n) (0) = 2n(2n—1) J^(2 Z) 2b ' 2 (?i 


21(271 3) (2 J g') 2rt ~ 4 tr 1 (2n_3) (^ 2 ) 


I (2n—2)(2n>—3)(2n—4)(2n—5) {21Q 2n “ 6 0^ 2n “ 4 *(if2) 


-L 


• (Ryshik and Gradstein 1957, 0-432) 
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where 


and so 


Gi{K*) = 


e K2 ‘J 0 (K2R) 

K*-K' 


G^ n) {K z ) = 27 (-1)' 

r = 0 


r! (JT*—JT , > -tr+1> 


~Tr \* 2 z Jo(K*R) 

d(K 2 ) L J 


n 


■* 2 (-l) f 

r=0 



r! (K2-K')- {r+1) R n - r J 0 {n - r \K*R) as 2 -» 0 - . 


Consequently, 


lim <? (2n) (0) = — 
0 - 


(2a)! V 
(»-l)l,.o 



lim J^"~ r ~ A \K 2 R), (45) 

*:->o+ 


on making /x 0 + . 

From (43), we find that dies out like t ~ 4 (or £“ 5 ) as £-> oo. We also note that 
when /x is not an infinitesimal, the law of decay of with time is an exponen¬ 
tial one. The use of /x thus physically signifies a quicker attainment of the 
steady state in presence of frictional forces. 


Analysis of or the steady state 
Lemma . 


J dK = - g [H 0 (Z^) + r 0 (Z' ff ) + 2tJ 0 m], 9 > 0, .. (46) 

where H 0 (x) is Struve's function of order zero. 

For 

J =f-r 0 (-*m ? > o 

(Erd&yi 1954, 14.3 (3)). 

Since 

H 0 (ze im7T ) = e^Hoiz) 

Y 0 {ze im7T ) = y 0 (z)+2iW 0 (z), m = 1, 2, 3,... 
the result (46) follows. 


(a) Concentrated pressure 
(i) Exact solution 

If 


Vo{x,y,t) =^~8(r)e < "', 

the wave-elevation £ 2 in the steady state is given by (17) 
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9 pU = ~ 2 ^ S ( r > e ‘ <a ‘~ 


K'P 

2n ( 


, f KJ 0 (Kr) 

Jo K ~ K ' 


dK 


P K'P K'ZP 

= - — 8(r)e-'-^ : e^+ — e-W*V)+ r 0 (^V)+2iJ 0 (^V)] 

by (46). 

On putting ft — 0, we have, for r > 0, 


(JJ^P . P(JL>^ 
» { «=-2^ e + 


^H*»(t)Wt) w (t)] " l4,) 


which gives the exact value of £ 2 for concentrated pressure, 
(ii) Asymptotic solution 


Cl) 

For — !> 1, the complete asymptotic development of £ 2 can be at once 

u 

obtained from (47). Retaining, however, only up to the second order terms, 
we have 


, _ p»*_(■_ jl\ -i(°r "“"■?) u.n lY^r 1 

gpU - (2nrg*)l \ + ° [U / . ' 


(48) 


The first order term of (48) agrees with that obtained by Lamb (1904, eqn. 
(135)) for this case. 

(b) Finitely extended symmetrical pressure . Exact solution 

If Po(*,y,t)=f(r)* iu *, 

where f(r) = 0 when r > a, 
the wave-elevation £ 2 in the steady state is given by (17) 

Kf(K)J 0 (Kr) 


K-K' 


where 


9PI% = f 

J o 

f(K) = f OLf(OL)J 0 (K*)dx. 

J 0 


dK , r > a 


On interchanging the orders of a- and if-integrations, we get 

gpU = —K’e iwt J a/(a)da J ^~Y f J o( K *) J o( Kr ) dK , r > a. 


Neumann's addition theorem for Bessel functions gives for r > a, 

/ o (^A)d0 


J 0 (Kx)J 0 (Kr) = - P 

o 

Jo(K«)Y 0 (Kr) = ± f 
Wo 


Y o (KX)d0 


(49) 

(50A) 

(50B) 
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where 


A 2 = r 2_j. a 2__2ra cos 6. ., 


• (51) 


Hence, after a change of orders of integrations in (49) by (50A), we get 
the A-integral of (49) = l - f ^ - y f [tf 0 (A'A)+ Y 0 (K'X)+2iJ 0 (K'X)]dd, (52) 

J 0 Jo 

on using (46) and the well-known result 


P° 

J o 


(KX)dK = j , A > 0. 


In (52), expanding H 0 (K' A) in a series of powers of K'X and integrating term 
by term, we obtain 


»TT QQ ^ 

J H 0 (K'X)d6= J (- 1 


> — C^ 2 -* 2 ) [r(»+t)]- 


-ft 


- «■ 2 <_1)n Rf( r2 -* 2 ) tnn+t)]-*P.+» 


r*+«2 2r« a \ n+i ja 

~9 -9 — -o cos 0) d0 

r 2 —a 2 r 2 —a 2 / 


- 2 d /r 2 +« 2 


r > a, 


.. (53) 

in which we have used Laplace’s integral for the Legendre function P y (z). 
This integral also gives 


CdQ it p /r 2 +a 2 \ 

J 0 A* (r 2 -« 2 )» -‘t 2 -* 2 /’ 


•• (54) 


Expressing the Legendre functions in (53) and (54) in terms of hypergeo¬ 
metric series in powers of ~, and evaluating the remaining integrals of (52) 

by means of (50A and 50B), we find 

the A'-integral of (49) = - 2 ^ (i J; 1,^) 


A' V , ,x„(AV/ 2) 2b+1 „ l _j 

2 Z, V mn+m z '- Fl \ n 2 ’ w ^ ’ 1 ’ r 2 / 


^ • [A»+l)] 2 " \ 

w ” U 

• ./ 0 (A'«)[ 7 0 (A'r)+2»V 0 (JTr)], r > a. 


Letting n -*■ 0 -j-, we get from (49) 


This gives the exact solution for £ 2 in the present case. 


„ (oj*rj2g) 


+£> — »+i; 1, 
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(c) Asymptotics for finitely extended pressure area 

If the impressed pressure p 0 (x , y , t) vanishes outside a finite region D 
of the xy- plane, we have, from (14), 

g P U = e Wt f j x, y outside D. 


On changing the order of the (K 1 , Af.>)-integral with that of the integral for F, 
we get 


,-W- a )+ K 2(y-m 


K-K' 


=-£>/>•» •• 


dK x dK 2 


(56) 


where 


R2 = (#—a) 2 + (i/—^) 2 . 

As in (5) above, 

1 ttK' 

the JUintegral of (56) = ~ ~ [H 0 (K'R) + Y 0 (K'R)+2iJ 0 (K f R)]. .. (57) 


In the asymptotic expansions that follow, we assume 



where 


r; = (« 2 +£ 2 ) max , r* = x 2 +y 2 . 


(58) 


Using now the asymptotic series for the different functions in the second 
term on the r.h.s. of (57), we obtain 


9pt 2 = 



T r F(ol, p)d<xdp+0 



Expanding R~*, R, R~* appropriately, we have, to the above degree of 
approximation, 


9pCz = 


———re *( e wt I) f f fl+i- {(a cos 0+|3 sin 0) 

(2irrg*)* J J D L 2r ( 

- ~ (-a sin 6+p cos 0) 2 +^ i |] F( a, fieT (acosfl+ ^ 3in 6> dadp 


- - (5)'"*'" & [?+A +° [(t) '] ,69) 


under the restrictions (58), 
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where 


When 

and D is a circle of radius a 
where 


s F (y cos-0, ~ sin 0^, F 99 = 


d-F 

W' 


F(x, y) = f(r) 

F(K 1 , K s ) = j(K) 

f(K) = f ccf(cc)J 0 (K<x) da. K = (Ai+Al) J . 

J 0 

Therefore, the above expression in this case reduces to 

When 

F(x, y) =fi(r) cos m0+f 2 (r) sin md r 
and D is a circle of radius a, 

F(K i3 K 2 ) = e im7T ’ 2 [Ji(iT) cos nufi+j 2 (K) sin rmf] 


where 


K=(K 2 1+ Kl)\ t = tan-1 


In this case, equation (59) gives 

y a/ 2xr V i ( U)2 ‘ r 77 mir\ 

x [?• (7) “ s ”»*+Js (y) si " ” 9 ] [ 1+ jis : < 1 

m. 


+0 


(d) The mean rate at which an impressed pressure 

COS 

Po(x, y, t) = F(x, y) (tu<) 

sin 

does work in generating waves can easily be obtained. 
From (14), we have 

gptz = -F(x, y)e«*~ ~ e«* f f F(«, 

" —00 

x ~ {H 0 (K’R) + Y 0 (K'R)+2iJ 0 (K'R )}] 

as in (57) above. Putting y. = 0, we get 
8 


• (60) 


(61) 


(62) 
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= ~[F(x, y)+G(x, 


F(z,P)J 0 ( c ^)<fodp .. (63) 


where 


° {x ’ y) = 2^ 




•(f)! 


a real expression. 


Taking the real or imaginary parts of ~ according as p 0 fF in (62) is cos cot 
or sin cot , we find, in either case, that the mean rate of work is 


= 


F (x> y) *’ I F(Z, P)J 0 I I dadfi 1 dxdy • 


'•(f)- 


-—— cos 

2 R\ If" f> a i a . ,0-jB 

— I = — I o dtp, 9 = tan -1 -—- 


it follows after a change of order of ^-integration, 


W ^ = Wp J (^ 1 +^ 2 ) #> 


where and i^ 2 are the real and imaginary parts of 


F j — cos tft, — sin iff 


F (x, y) = /(r), 


W n „ = 


7T 2 0> 5 ‘ /CU 2 \] 2 


which is the expression derived by Wehausen and Laitone (1960, 21.9). 


Particular cases . Central pressures over elliptic regions 
Some of these are defined by 


where 


Po(*> V> 0 = F (*> y) M 

sin 


x 2 y 2 \ v - 1 


F{x >y )=S Mb V-a 2 ~b 2 ) ’ in8ide a 2+ b 2=1 


x 2 y 2 


, outside ,, 


For v = l,f, 2, these disturbances may be characterized as (elliptic) cylindrical, 
8b 
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ellipsoidal and (elliptic) paraboloidal pressures respectively, P being the total 
pressure on the surface in each case. We have 


F(K U K«) 


(A) Unsteady waves 


(.Ky+K 2 2 b 2 y 12 


( 66 ) 


(Erdelyi 1954, 1.3 ( 8 ), 1.13). 


For a pressure F(x, y) cos ojt, we have, from (39), 

[ J ^* E) cos (rY2 >~ 2 7F2 ain {/t ' 2) 


x j (t' 2 E) 2 ( et 2e -2.) J v+ - 2 (t'*E) 

+t' 2 E ^+ eZ ^^y y+ ,(t' 2 E) + hh(t'*E ){] 


under the restrictions (32), 
where 


t' = = (1 —e 2 sin 2 6)\ 

2 r a 


(67) 


.. ( 68 ) 


e being the eccentricity of the ellipse. 

For a pressure F(x y y) sin c ot, we have, from (39), 


np(ag) i r'(t' 2 E) 1 ' 


J v (t’ 2 E) sin (r'<' 2 )+ 2^72 cos 


(r'r-) 


X { (t'*E) 2 -2) J v+2 (t' 2 E) 

+t' 2 E (.5 + e2 - g- 2ff ) J v+ ,(t' 2 E) + lJ v (t' 2 E)\] .. .. (69) 

under the restrictions (32). 

When a = b, we get the corresponding expressions for the symmetrical pres¬ 
sures over circular regions. 


(B) Steady waves: 

(i) Exact solution for symmetrical pressures . 
Let 


Therefore 


Po(z, V: 0 = 5 I 1 -5) e ia “,r<a 


/(*) = 


= 0 , r > a 

2*’ 1 .'P1» J„(Ka) 


• (Kay 


( 70 ) 
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Also, 

J «/(*'2-F’i (—«+i —w+i; = -2-^1 (~ w +i - w +i; l ’+ 1 > " 2 ) 

(Erdelyi 1953, 2.4 (2)). 

Consequently, the exact solution for the wave elevation in the steady state is 
given by (55) 

=-1 fet (fP : 

X 2 F 1 (-«+i -« + *; v+1, j)+2‘'-2a-'vPr(v)(y) 

x Mt1 +w »(t)Mt)' , ' > ” .mi 


(ii) solutions for central pressures on elliptic regions . 

For a pressure F(a;, cos <»t, F{x , ?/) being as defined in (65), we have, from 

m, 

gp{<i = — - cos (rt\— <*>t+^\ -\ --—- sin ir'f — a>t+j\ 

H (irar'grfifaYl \ 1 4 / 8 r’EH\ \ 1 4 / 

x j (el,) 4 sin 2 20J v+ 2 (Et^+lEiettf cos 20J„ +1 (Et^)+E^J v (Et\) Jj .. (72) 


under the restrictions (58), 
where 


2 _ 

1- T~ 


(73) 


r', e and E being the same as before. 
For a pressure F(x , ?/) sin we get 

0 v-i 


9 Pt 2 =7— , I 'v!/p V sin co<+j)--—2 cos ( r t \~ wt +~i\ 

2 (T<ir’gZf(Ety [ \ 1 4 / Sr '.® 2 * 2 \ 1 4 / 

X {(e<i ) 4 sin 2 20«/ ( , +2 (£'<f)+ 4 £(c < 1 ) 2 cos 20 J ,. +1 (Et\)+E*J v (E%)] j .. ( 74 ) 

under the restrictions (58). 

When a = b, we get the corresponding expressions for the symmetrical 
pressures over circular regions. 


Discussions 

The similarity of the real and imaginary parts of (39) with the expressions 
for the elevation of Cauchy-Poisson waves due to an initial surface elevation 
and impulse suggests that the unsteady waves given by (59) behave, in many 
respects, like the corresponding Cauchy-Poisson waves for large values of r and 
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t and for certain frequencies determined by (32). This similarity is further 
evident in the law of decay of the unsteady state as shown in equation (42). 

As to the effect of the frequency, we observe (eqn. (39)) that, in the range 
of approximations (32), the amplitude of the unsteady waves varies as the 
frequency or as its square according as the impressed pressure has the phase 
sin cot or cos ot. 

At large distances from the seat of the disturbance in the steady state, the 
third and fourth terms of the expansion (60) for a symmetrical pressure 
f(r) e lix>t over a circle of radius a are respectively 


9 i 

l28co 




and 


_ o i0>l f a. 

w 2 r3 J Q 


if (cl) dec . 


(75) 


If, now, o and a be such that J j = 0, equations (60) and (75) show that 
the fourth order term is dominant and this term represents a standing wave of 
amplitude O^j as r-> oo. Further, equation (59) shows that this phenom¬ 
enon does not generally occur in motions due to asymmetrical pressure 
distributions. 

Lastly, we note that in any type of pressure distributions, as equation (59) 
shows, the amplitude of the diverging progressive steady waves of length 

dies out like as oo. 
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